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Abstract 

A precise study within variational approach of the basic properties of the three-particle spectrum 
and structure functions with Gaussian potential near the critical coupling constant of interaction 
where the Efimov effect takes place is carried out. A method is developed to calculate highly excited 
states with very small binding energies, and numerical analysis is carried out for the ground and 
three excited energy states. For these states, one-particle density distributions, formfactors, pair 
correlation functions, and momentum distributions are calculated. It is found that the second 
excited state has already all the basic features of a level from the infinite Efimov series. An 
essential asymmetry is found in the position of energy levels with respect to the critical constant 
point. A halo-type structure is revealed in the one-particle density distributions, and formfactors 
are shown to have specific dips of finite depth, the number of dips being equal to the number of 
the state. The behaviour of pair correlation functions and momentum distributions is studied for 
three-particle states. 
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I. INTRODUCTION 

As known, the Efimov effect Ql reveals itself in a system of three particles with finite- 
range interaction as the appearance of the infinite number of levels near the threshold with 
very small energies under the condition that the two-particle subsystems have an infinitely 
small binding energy or an infinitely large scattering length. In the case of three identical 
particles, the energy ratio for neighboring highly excited near-threshold Efimov states 0, 0] 
is A^; = lim (En/En+i) ~ 515. In three- nucleon systems, the conditions necessary for the 

n— >oo 

effect to take place can be easily broken by taking into account the spin structure of the 
nuclear interaction P, |^ or the Coulomb long-range repulsive potential. 

Theoretical studies of the Efimov near-threshold energy spectrum were carried out both 
by asymptotic expansions with the use of the small parameter rn/a <^ 1 (being the ratio 
of the range of forces to the two-particle scattering length) l|, Q, Q, S, Isl and by numerical 
calculations [3, 0] based on the Faddeev equations by using separable interaction potentials. 
The main properties of the energy states are universal to the great extent and do not depend 
on a specific form of the potential. For the special choice of an interaction potential in the 
form of two components with essentially different radii, new additional excited energy levels 
appear in the three-particle system — the "trap" effect In the case of local potentials, 
numerical calculations of the Efimov energy spectrum were not carried out systematically 
because of essential difficulties due to the different scales of very small near-threshold excited 
states energies and the ground-state energy. This requires to develop the precise methods of 
calculation of the highly excited weakly bounded states of three-particle systems with local 
interactions. 

Due to the use of Gaussian bases in the variational calculations of a few-particle sys- 
tems Q|, noticeable progress was achieved in the calculations of the main properties of 
few-body systems Inl H y Q with interactions of different nature in recent years. In 
the present paper, within the precise variational approach, a study of the near-threshold 
Efimov states of a three-particle system with Gaussian potential is carried out up to the 
third excited level, and all the main structure functions for various energy states are found. 
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II. STATEMENT OF THE PROBLEM 

Consider a system of three identical particles 

i=l i>j=l 

with a two-particle interaction taken in the Gaussian form for simplicity, where Vq is the 
intensity and tq is the radius of forces. Let us study the properties of the energy spectrum 
and the main structure functions of the system in symmetric (with respect to permutations 
of particles) states with zero total angular momentum in the range of parameters where the 
Efimov effect takes place. It is convenient to use the dimensionless variables and to measure 
all the distances in tq and energies in ff/mrQ. In this case, the dimensionless Hamiltonian 

3 3 

^ = -^E^^-^E-p(-4) (2) 

i=l i>j=l 

contains only one combination of physical parameters, the coupling constant g = mrQVo/ff, 
which determines the properties of the spectrum and wave functions of the system. 

The solution of the problem on eigenvalues and eigenfunctions for bound states is carried 
out in the framework of the variational Galerkin method in the Gaussian representation with 
the use of special optimization schemes (see [3,0, Q|) which enhance the convergence and 
ensure a high accuracy of calculations with the least dimensions of variational bases. Within 
this approach, the variational wave function of symmetric states of the three-particle system 
with zero total angular momentum can be presented as the expansion in Gaussian functions 
depending on relative distances = |rj — rj|, 

^ (ri2, ri3, r23) = ^ -DfcS ) = ^ D^S exp {-aurl^ - htrl^ - Curl^) , (3) 
fe=i fc=i 

with a set of variational parameters a^, hk-, c^. For the bound states, the linear coefficients 
Dk are found in the variational approach from the system of linear algebraic equations of 
a problem on eigenvalues (which is equivalent to the Schrodinger equation in the chosen 
Gaussian basis): 

K 

^Z}„((S^fc|iJ|S^„)-E(S^fc|S^„)) =0, fc = l,2,...,/C. (4) 

n=l 
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Note that all the matrix elements in Eq. (0} for Hamiltonian Q can be calculated in explicit 
form. The chosen method of calculation has proved its high efficiency and accuracy in a 
number of problems of a few particles with interactions of different nature. But, in the case 
of the near-threshold Efimov spectrum, we are faced with an additional difficulty due to very 
small values of the energies of these levels. Since the ratio of the Efimov neighboring energy 
levels \e = lini (En/En+i) ~ 515 is rather large, there exists a real difficulty for all the 

n— >oo n 

methods (including the method based on the Faddeev equations j9| and variational methods 
as well) in calculating the energies with essentially different scales and the corresponding 
wave functions with oscillations at essentially different distances. 

In the present work, the main structure properties of the ground and three excited states 
of a system of three particles are calculated with the coupling constants near g ^ Qcr, where a 
two-particle bound state appears and the Efimov effect reveals itself. The numerical analysis 
of highly excited near-threshold Efimov levels appears to be a nontrivial problem for all the 
methods including the variational method because of the necessity to prepare the "initial" 
configurations for the wave functions of highly excited near-threshold bound states to start 
the minimization of energy in the nonlinear variational parameters , bk, Ck- 



III. THE PROBLEM OF INITIAL CONFIGURATIONS 

An advantage of the variational method lies in the possibility to calculate (in the approx- 
imation of a given number of the basis functions) the whole spectrum of low-energy bound 
states for a given Hamiltonian or to calculate the energy and wave function of the chosen 
ground or excited state separately. The further minimization in the nonlinear parameters 
(^k, bk, Ck enables one to improve the variational estimation for the energy and to achieve the 
best result at a fixed dimension of the basis in (jS)). It is clear that if the initial configuration 
gives, in fact, the bound (ground or excited) state, we have to approach the exact solution 
for this state by minimizing the nonlinear parameters a^, bk, Ck and expanding the basis 
dimension. However, at a finite and even rather large number K, of the basis functions in Q 
and with a random set of the nonlinear parameters a^, bk, Ck (we call it a random "initial 
configuration"), nobody can be sure not only in the fact that the studied energy level is 
close to the exact value, but, even more, that it is below the threshold of a decomposition 
into subsystems at all. In the latter rule, the further minimization in nonlinear 
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parameters should be terminated on approaching the rated level to the threshold (from the 
side of the continuum spectrum), because the wave function can have nothing in common 
with a bound state. Such problems reveal themselves starting from the second excited state 
which exists in a narrow interval near g g^^. = 2.684005 (the critical constant for the 
chosen Gaussian potential). There is no such problem for the ground state of three parti- 
cles, whose initial configuration can be simply prepared even with one Gaussian function at 
sufficienly large g. The mentioned problem is also easily solved for the first excited state 
if one takes several Gaussian functions, among which at least two or three ones depend on 
the parameters ak-,hj^ Ck taken from the ground-state function, while the rest have to form 
a cluster function [l3[. This enables one to rather simply construct the necessary config- 
uration (even in the case where the initial configuration corresponds to an energy slightly 
exceeding the two-particle threshold). For the ground and first excited states, the problem 
is simplified also due to the fact that these levels exist at all the interaction constants g 
greater than certain critical values (see Tabled)). Therefore, it is sufficient to enhance the 
attractive potential, to form the initial configurations at a sufficiently large g, and then, by 
gradually decreasing g (the evolution in coupling constant), to achieve the necessary value 
of g by minimizing the energy in the parameters a^, bk, Ck- We used both variants to form 
the initial configurations for the ground and first excited states. 

For the following excited near-threshold Efimov states with extremely small binding en- 
ergies existing within a narrow interval of g near the critical two-particle constant gcr, the 
problem of initial configurations requires special methods for its solution. One way is to 
use a wider class of Hamiltonians, in particular, those with an additional interaction giving 
a possibility to form relatively easily the configurations of highly excited levels. Then the 
above interaction is eliminated step by step (in order to return to the initial problem) with a 
simultaneous change in the parameters a^, bk, Ck within the variational procedure in order to 
keep the energies of the formed configurations below the threshold. The Coulomb attractive 
interaction seems to be suitable for this purpose giving rise to the infinite number of levels 
for the three-particle system. An additional advantage of this interaction is the fact that all 
matrix elements in the Gaussian basis for the Coulomb potential are known to be calculated 
in explicit form. 

Another way to prepare the initial configurations is to deal with a Hamiltonian with 
different masses of particles. It is known that the Efimov effect takes place in three-particle 
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systems with different masses Moreover, the ratio of the energies of neighboring 

levels depends significantly on the ratio of masses (see also ^^). From the analysis of the 
power asymptotics of a three-particle wave function in the momentum representation within 
the approach based on the Faddeev equations, it can be shown that the following secular 
equation is the condition for the solvability of three-particle equations: 

L(l; 2, 3)L(2; 1, 3) + L(l; 3, 2)L(3; 1, 2) + L(2; 3, 1)L(3; 2, 1) + 



L(l; 3, 2)L(2; 1, 3)L(3; 2, 1) + L(l; 2, 3)L(2; 3, 1)L(3; 1, 2) = 1; (5) 

where 

■ / ^ ATfl^i'' Aijfcx r^ sinh(so6'/2) a - 6^ 

L{i;j,k) = N{^;^); N{a]h) = -—= — — ; cos^ = — — ; 6 

IXi ruk Oi/a So cosh(sovr/2) a + 

= ; /i* = — ^— . 7 

rrij + rUk mi + rrij + 



In the case of equal masses, the secular equation (jH)) becomes essentially simpler 



iV(^;i)4. (8) 



or 

8 sinh(7rso/6) 



1 (9) 

SoCosh(7rso/2) 

and one has for the power sq of the asymptotics of the Faddeev equations solution (the 
Danilov-Minlos- Faddeev- Efimov constant), Sq = 1.0062378. This constant determines, in 
particular, the ratio of Efimov neighboring energy levels (the scale invariance multiplier), 

\e = lim {En/En+i) = exp ( — ) ^ 515.035 . (10) 



Vso 

Fig. d shows the dependence of the ratio A^; = exp(27r/so) of the energies of Efimov neigh- 
boring levels on the ratio of the masses of particles. It is seen clearly that the case of equal 
mases is the most unsuitable for the Efimov effect to reveal itself (the scale multiplier A^; is 
the greatest) and, at the same time, it is the most difficult case for numerical calculations. 
But, in the case of essentially different masses, one may have the ratio of energies not too far 
from 1 instead of (10) (for any ratio of masses, of course, we have A^; > 1 ). In particular. 
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in the hmit mi/m^ = m2/m^ (the one-center problem), relation yields the limiting 
secular equation for sq: 

{l + 8 [l + cosh-^ i^-f)] [sinh (^) /so] '} cosh- (^) = 1 • (H) 

This results in so = 1.139759, and the corresponding ratio of energies A^; = 247.826. If, vice 
versa, mi/m^ = 1712/1713 —>■ 00 (the two-center problem), the value of sq increases infinitely 



as So — ^ k^ymi/m^, where the coefficient k = 0.40103 is found from the secular equation 

exp(-W2)_ 

kV2 ~ ^ ^ 

The corresponding ratio of energy levels goes to 1 as Aij ^ 1 + at ma/mi 0. It 

should be noted that, due to the essential dependence of the parameter A^ on the ratio of 
masses, it is more probable to observe the Efimov effect in a real physical system just with 
different masses of particles whose pairwise interaction is of the resonance type. Close to 
such systems might be three-cluster nuclei consisting of a rather heavy magic core and two 
neutrons. In particular, if the nucleus ^°Ca serves to be the core, it follows from (0) that 
A^; ~ 262 (i.e., it is twice smaller than that for equal masses). If, to say, we consider the 
system of mesoatoms of deuterium and tritium interacting with an electron, then, due to 
the small mass of an electron in comparison with the masses of mesoatoms, one has a large 
value of So from (0) enough to result in A^; = 1.26. This value is not essentially greater 
than 1 and much less than that in the case of equal masses. This means that the energies 
of neighboring Efimov levels would be of the same order of magnitude if they were observed 
in such a system. 

Returning to the problem of formation of initial configurations, we note that the configu- 
ration for a given (the n-th excited) state can be prepared rather easily at different masses. 
Then one has to restore step by step the equality of masses by minimization of the energy 
of this level in nonlinear parameters. In this way, we can get the necessary variational wave 
function corresponding to the configuration of a bound state. 

To obtain the initial configurations for variational wave functions, we paid the main 
attention to the method of evolution in the coupling constant g at a fixed (zero) energy, 
which is possible because the interaction potential in ^ has a definite sign. The method, 
in essence, consists in a variational procedure with respect to the critical coupling constant 
g to be determined for the given ground or excited state at a fixed (zero) energy. It is 



the general statement that the variational method can be formulated with respect to the 
coupling constant g at a fixed energy provided that the interaction potential is negative 
definite. Taking an arbitrary initial configuration with a sufficient number of the basis 
functions with some parameters a^, bk, c^, we find the corresponding coupling constant 
g of the chosen excited state by variational procedure and advance step by step towards 
small g by varying the parameters a^, bk, Ck (or sometimes, if necessary, by spreading the 
basis). Thus, we approach the two-particle critical constant from the upper side, where a 
two-particle ground state appears. When we are close enough to the region g ^ gcr, where 
the Efimov excited state under consideration must exist, we return to the commonly used 
variational procedure in energy using the obtained set of the parameters a^, bk, and start 
to increase the parameter g. If we were close enough to the critical constant gcr, the three- 
particle energy level crosses the two-particle threshold and goes down below the threshold. 
Thus, we obtain the required initial configuration. Such a possibility is connected with the 
fact that the two-particle energy threshold behaves itself by the square law near the critical 
coupling constant, 

E(2) = - 0.107 (13) 
while three-particle levels reveal a linear behavior Q|, 

En = Ki9- gtjt) , (14) 

where is the critical constant for the n-th level appearance on the left from the critical 
constant, and kn = dEn/dg (at the critical point g = g[^j^) is the corresponding angular 
coefficient (see Table H}. Using the above method, we succeeded in preparing the initial 
configurations for the second and third excited states which exist near the critical constant 
gcr and belong to the Efimov series of levels. Note that the greater the number of the 
excited near-threshold Efimov level, the greater are the difficulties in preparing the initial 
configuration. It is not only due to the necessity to increase the dimension of the variational 
basis, but also because the nonlinear variational parameters differ from one another by 
several orders of magnitude. Moreover, there appears a considerable quantity of nonlinear 
parameters such that they are roughly by times smaller than those typical of the previous 
energy level. This makes it necessary to significantly decrease the steps in the process of 
variation in nonlinear parameters and decelerates the whole procedure of preparing the initial 
configuration. Because the energy of already the third excited level is less by nine orders 
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of magnitude in comparison with the ground state energy and the dimension of the basis 
which is necessary to ensure a rehable accuracy achieves several hundreds of functions, the 
calculation of levels higher than the third one becomes a serious problem. 



IV. ASYMMETRY OF THE EFIMOV ENERGY LEVELS 

It is convenient to depict the calculated energies of the ground and excited states of 
a three-particle system as the square roots of the moduli of three-particle energies minus 
the two-particle energy threshold versus the reciprocal scattering length (mstead of the 
coupling constant g), because such dependences are model-independent lilslQ, to the 
great extent. Since the orders of magnitude of the values (of both energies and reciprocal 
scattering lengths) are essentially different for different levels (the ratio for neighboring levels 
being about \/Xe = exp(7r/so) = 22.692 ), Fig. |21 shows each level in its own scale in order 
to emphasize the universal properties of highly excited states. Namely, the horizontal axis 
depicts a~^A^~^^''^, and the vertical axis shows (i?(2) — -E„)^/^A^~^''^^. The ground and 
first excited states exist at any g greater than the corresponding critical values (see Table 
H]). The reciprocal scattering length is approximately proportional to g — gcr near the 
two-particle critical constant g = gcr- Note that the separation of three-particle levels from 
the zero energy occurs by law (jl4p . though this linear behavior takes place only in a very 
narrow interval near the corresponding three-particle critical coupling constant, and the 
angular coefficient kn decreases essentially with increase in the number of the excited state. 
The second and third excited states appear with increase in g (in the region g < g^r) and 
then disappear (in the region g > gcr) at the two-particle threshold (the corresponding values 
of the coupling constants g and the reciprocal scattering lengths are given in the second and 
third columns in Table|l}. A rather good coincidence of the curves for the second (n = 2) and 
third (n = 3) excited levels depicted in Fig. El testifies to that the asymptotic formulae 
are valid already starting from the second excited state and indicates the fact that highly 
excited Efimov states possess the scaling property. Although we restricted ourselves by the 
calculations of only three excited levels, it should be assumed that the rest of levels is to be 
determined by the asymptotic relations, and thus the properties of the whole spectrum are 
completely described. A significant asymmetry of the curves obtained in our calculations 
for the second and third excited states relative to the critical concentration point gcr of the 
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Efimov spectrum testifies to the asymmetry of all the rest levels. Each of them extends to 
the right side from the critical point much farther than to the left. This essential asymmetry 
of the levels relative to the point g = Qcr (or a^^ = 0) is connected with the following fact. 



ranges up to 



Though the potential of the effective long-range interaction ~ —R~^ 
distances of the order of \a\, it remains attractive at essentially larger distances at a > 0. 
Whereas, at a < 0, it becomes practically zero out of the region of about \a\ (see also jisl]). 
The above-mentioned asymmetry is clearly seen in Fig. El where the second and third excited 
states are depicted in natural scale versus the coupling constant g. This figure contains a 
fragment of the near-threshold area in larger scale in order to make the third level visible, 
since it looks almost like a dot and cannot be distinguished from the two-particle threshold 
within the main figure. Note that the largest binding energy of each of the Efimov levels 
(minus the energy of the two-particle threshold) happens at the constants rather far from 
the critical concentration point of the spectrum (see Fig. |21and Fig. EI), though the infinite 
number of levels appears only at the limiting point Qcr- The three-particle levels cross the 
two-particle threshold (with increase in g) at very small angles at the points in Fig. 01 where 
the solid lines of these levels become the dashed line of the two-particle threshold. After 
crossing the two-particle threshold, the three-particle levels may exist over the threshold 
as virtual states (see It is worth to note that, with the essential enlargement of the 

(2) 

coupling constant g up to g > gcr = 16.3, the second excited level appears below the two- 
particle threshold, the third one does a.t g > gc? = 30.8, the fourth does a.t g > gir = 32.55, 

(5) 

the fifth does at (? > gcr = 50.89, and so on. Their binding energies increase with the 
coupling constant. The angular coefficients determining the angle between the curve of the 
n-th level and the two-particle energy threshold (near the critical point of the corresponding 
level appearance), 

E^-E(^2) = ~kn{g-gi'^^) , (15) 

are as follows: h = -0.0045, h = -0.1, h = -0.01, and h = -0.03. The fact that, for 
the third and fourth levels, gif^ and gif^ are close one to another, as well as that the angular 
coefficients and ^4 are of different order of magnitude, testifies to a different nature of 
these levels. 
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V. STRUCTURE FUNCTIONS OF THE EFIMOV STATES 



The wave functions of the three-particle ground and excited states calculated within 
the precise variational approach enable us to calculate directly such structure functions 
of the system as density distributions, formfactors, pair correlation functions, momentum 
distributions, etc. It is a rather simple problem to calculate any such an average since 
the wave functions of each state have the simple form (jS)) of a superposition of Gaussian 
functions with the already known linear (Dk) and nonlinear (a^, bk, Ck) parameters. 

Fig. m presents the one-particle density distributions, 

^W = (^^|:5(r-(r,-R_.))^ , (16) 

for the ground and excited states at the critical two-particle coupling constant g^r- The 
density distributions (fT^ are shown in logarithmic scale on both axes. This is made, on 
the one hand, in order that the typical "halo" structure appearing for different states 
at distances of different scales can be shown in the same figure. On the other hand, 
on the chosen logarithmic scale, the density distributions reveal clearly the almost peri- 
odic dependence of the wave functions of Efimov states on the global radius logarithm ^ 
~ sin ^|so| In {^^Jra \E\ /h?^ + A^. Moreover, this asymptotic behavior is valid from the 
distances of about the radius of forces (in our case, of about 1) to the distances of about 
h{m\E\) ' . The less the binding energy of an excited state, the longer is the extension 
of this asymptotic behavior. The density distribution of the first excited state is seen to 
change sharply its behavior between the short and long distances. The long-range "halo" is 
situated around the more condensed central core similar to the density distribution of the 
ground state. The density distribution of the second excited state changes its behavior two 
times having a "halo" in the form of two concentric spherical layers put one into another. 
They are of essentially different radii and densities and have the central part of the same 
type as in the case of the first excited state. The density distribution of the third excited 
state changes its behavior three times (three concentric layers around the central core), and 
so on. The ratio of the radii of neighboring layers is about ^/Xe — 22.7. Since the density 
distributions are normalized by J n{r)dr = 1, they decrease at short distances when the 
number of the excited state increases, because the size of the system is growing. 

The fact that the ratio of the sizes of neighboring states of the system is about \/Xe is a 
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consequence of the universal model-independent properties of Efimov states and is confirmed 
once more by calculations of the r.m.s. radii (see Table U with the calculated Rrms at the 
critical value Qcr)- Note that the smallest value of Rrms for a given Efimov level is achieved 
at the constant g corresponding to the largest binding energy, rather far to the right from 
the critical point gcr (see Fig. |2I). When the coupling constant approaches the value, where 
the Efimov levels appear at the threshold (i? — > 0), the calculated mean square radius 

2\l/2 ^ ^ 



Rrms ^ = rMr)drj (17) 

reveals the behavior 

(r^>^/^ = -^ + c, + .... (18) 



That is, the radius goes to infinity, in accordance with the general physical principles, as 
(r2)^/2 ~ (E)-i/2 at E ^ (we recall that the energy level near the threshold depends 
linearly on the coupling constant according to (jl4j) ). It is important to notice that the 
growth of the r.m.s. radii at ^ is connected mainly with the extension of the region 
where the asymptotic behavior of the wave function ~ sin ^|so| In (^\Jm \E\ /h?^ + is 
valid, i.e., up to the distances of order h{m \E\) . At the same time, the wave function 
at short and intermediate distances does not change (except for the normalization factor 
depending on the total size of the system). As a result, in all the intervals of g where a 
given Efimov level exists, the behavior of the wave function and, in particular, the density 
distribution is not changed at short and intermediate distances (but only the normalization 
does), which is confirmed by Fig. |S1 This figure shows the density distributions of the second 
excited state at several different values of the coupling constant, from the point g = 2.682195 
near the level appearance to the point gl./gf^f = 2.7244 of its disappearance. Other typical 
constants chosen by us are as follows: the two-particle critical constant gcr (the point of the 
Efimov spectrum concentration) and the constant g = 2.702, where the energy of the second 
excited state (minus the two-particle threshold energy) is the greatest. 

Specific sharp changes in the density profiles near the distances, where the corresponding 
wave functions pass through zero, manifest themselves in formfactors, 

-^(9) = J exp(— zq • r)n(r)(ir , (19) 

in the form of specific dips of finite depth. Fig. IHl shows the formfactors of the ground state 
(no dips), the first (one dip), second (two dips), and third (three dips) excited states. A 
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logarithmic scale commonly used for formfactors on the vertical axis is accomplished with 
the logarithmic scale in the squared momentum transfer in order to show clearly the almost 
periodic repetition of the dips connected with regularities in the "halo" structure of the 
density distributions of excited states. Since the greater distances in a density distribution 
correspond to a less momentum transfer in the corresponding formfactor and vice versa, a 
new dip in the formfactor of the next excited state appears at a less momentum transfer, 
whereas the previous dips remain practically at their positions. 
Fig. [7| shows the pair correlation functions 



of the ground state and three excited states. We multiply the pair correlation functions by 
in order to show clearly the following interesting effect. With increase in the number of 
the excited state, we observe the spreading of the area, where the pair correlation function 
is proportional to the reciprocal squared distance between a pair of particles. That is, the 
two-particle subsystem in the system of three particles has the density distribution very 
similar to the squared two-particle wave function. This gives rise to the effective long-range 
interaction ~ — R in the system of three particles in the region from distances of 

order of the radius of forces i? ~ tq to the distances of about R h{m where the 

wave function starts its exponential decrease due to a finite binding energy. 
The momentum distributions 



for the ground and three excited states are presented in Fig. |H1 It is found that the mo- 
mentum distributions vanish at approximately the same limiting momentum po (in our case, 
10° < Po < 10^). This is due to the fact that the largest momenta in the considered three- 
particle system are typical of the shortest distances in the coordinate representation (of 
order of the radius of forces). There, the wave function has the sharpest behavior, but it 
is practically unchangeable under growing the number of an excited state. With increase 
in the number of a state, the wave function extends to larger and larger distances almost 
without change in the central part (see, for illustration. Fig. |3] depicting the density distri- 
butions). Note that, with increase in the number n of the state, the momentum distribution 




(20) 




(21) 
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n 



D{p) increases essentially at small momenta proportionally to |^A^ j . Indeed, since the 
size of the system increases by \/\e times for each of the next excited states and its vol- 
ume increases by (A^;)^^^ times, the region of momenta where the momentum distribution is 



get that the momentum distribution density grows by {XeY'^ ^ 1.17 X 10^ t imes at small 
p with every increment of the state number. Fig. |H1 also demonstrates the extension (to- 
wards zero momentum) of the area of the power dependence of the momentum distribution 
D{p) ~ which looks like a linear decreasing function on the logarithmic scale. Such a 
behavior is connected with the fact that the greater the number n of the state, the greater 
is the role of the kinetic energy entering the Green's function in the formation of the power 
dependence of the wave function in the momentum representation. Really, in the formal 
operator representation of the Schrodinger equation solution 



the energy ^ with increase in n. In this case, instead of the Green's function, we have 
only the kinetic energy in the denominator. The expression V \^n) varies slightly, because 
it is nonzero only within the region of the short-range interaction, where the wave function 
is changed insignificantly at growing n. After the integration with the 6 - function, the 
squared modulus of the wave function in D{p) leads to the value in the denominator which 
is proportional to the squared kinetic energy of one particle, by giving rise to D{p) ~ 

VI. CONCLUSIONS 

We have developed a variational approach with the use of the Gaussian basis and precise 
optimization schemes and have studied the spectrum and wave functions of three particles 
with Gaussian pairwise interaction near the critical coupling constant, where the Efimov 
effect takes place. A special technique is developed for preparing the initial configurations 
of highly excited near-threshold weakly bounded energy states using both the deformation 
of the initial Hamiltonian (by changing the interaction or the mass ratio) and the method 
of evolution in the coupling constant at a fixed (in particular, zero) energy. For several 
calculated levels from the Efimov series, it is shown that they have an essential asjTiimetry 
with respect to the critical point and the largest binding energy (minus the threshold energy) 





(22) 
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significantly far to the right from the critical coupling constant which is the point of the 
spectrum concentration. The scale invariance is confirmed for the Efimov levels starting 
from the second excited state of three particles. 

For the first time, we have carried out the study of the specific behavior of the density 
distributions, formfactors, pair correlation functions, and momentum distributions for the 
ground and three excited states of the Efimov series. The details of the density distributions 
are found, and the halo-type structure is revealed for the excited states. It is shown that the 
corresponding formfactors have some dips of finite depth which are arranged periodically in 
the momentum transfer on the logarithmic scale. The behavior of pair correlation functions 
and momentum distribution profiles is found, and the main conclusions following from the 
asymptotic formulae for the Efimov spectrum are confirmed. The main specific properties 
of the structure functions of three particles within the area of parameters where the Efimov 
effect reveals itself would take place with other attractive interaction potentials as well. 

The precise calculations of the spectrum and the wave functions of highly excited near- 
threshold levels confirm the high accuracy of the variational method in the Gaussian repre- 
sentation combined with the efficient technique for the maximum optimization of the basis. 
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TABLE I: Characteristics of three-particle energy levels in the Efimov effect region. In the second 
and fourth cohimns, the corresponding values of the reciprocal scattering length are given for 
each value of the coupling constant g. 



n 


(n) 


kn = dEn/dg 


(n) 
dright 

(a-M 


En 

(a-1 = 0) 


(a-1 = 0) 





2.13079 

(-0.22853) 


-0.147 




-0.23845 


1.21 


1 


2.64396 
(-1.35046x10-2) 


-4.72x10-2 




-4.514x10-^ 


2.25x10^ 


2 


2.68215 

(-6.1717x10"'^) 


-2.07x10-^ 


2.7244 

(1.32485x10-2) 


-8.76x10-'^ 


5.04x10^ 


3 


2.68393 
(-2.48262x10-5) 


-1.38x10-5 


2.68567 
(5.53501x10-^) 


-1.7x10-9 


1.14x10^ 
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FIG. 1: Coefficient A^; for the neighboring Efimov energy levels versus the mass ratio (ei2 = mi/m2, 
and £23 = 1^2/1713) for the system of three particles. 
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FIG. 2: Dependence of (£^(2) —En\^l'^ on the reciprocal scattering length oT^ for the ground (n = 0) 
and excited (n = 1, 2, 3) states. Each of the levels is depicted on another scale with the use of the 
multiplier A^"^''^^ for both axes. 
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FIG. 3: Dependence of the second and third excited energy levels on the coupling constant g. The 
third excited level is shown in the enlarged fragment. 
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FIG. 4: One-particle density distributions for the ground (n = 0) and excited (n = 1,2,3) states. 



21 



n(r) 
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FIG. 5: Density distribution profiles with several different coupling constants g for the second 
excited state. 
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FIG. 6: Formfactors for the ground (n = 0) and three excited (n = 1, 2, 3) states. 
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FIG. 8: Momentum distributions for the ground (n = 0) and three excited (n = 1,2,3) states. 



25 



